Non— Interacting Fermions in a One— Dimensional Harmonic Atom Trap: Exact 
One— Particle Properties at Zero Temperature 
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One-particle properties ol non-interacting Fermions in a one-dimensional harmonic trap and at 
zero temperature are studied. Exact expressions and asymptotic results for large Fermion number A*' 
are given for the particle density distribution no^z, N). For large A*' and near the classical boundary 
at the Fermi energy the density displays increasing fluctuations. A simple scaling of these tails 
of the density distribution with respect to N is established. The Fourier transform of the density 
distribution is calculated exactly. It displays a small but characteristic hump near 2kF with kp 
being a properly defined Fermi wave number. This is due to Friedel oscillations which are identified 
and discussed. These quantum effects are missing in the semi-classical approximation. Momentum 
distributions are also evaluated and discussed. As an example of a time-dependent one-particle 
problem we calculate exactly the evolution of the particle density when the trap is suddenly switched 
off and find a simple scaling behaviour in agreement with recent general results. 



PACS Nos. 05.30FK, 71.10.Pm, 0375.Fi 
I. INTRODUCTION 



Recent years brought about spectacular successes in 
the study of dilute bosonic quantum gases confined to 
atomic traps at extremely low temperatures. These and 
the experimental details are reviewed in ^. 

The next stage of investigations will incorporate 
fermionic quantum gases. Fermi degeneracy of potassium 
atoms ("^^K) has recently been observed in Q. The effects 
of interactions between the neutral atoms are of particu- 
lar interest. They can give rise to collective ground states 
like superfluid phases. 

Another development regards the construction of 
highly anisotropic traps, e.g., the microtraps in P-u). 
The magnetic trapping fields can be taylored so as to 
make the confining potential harmonic. 

If the longitudinal confinement frequency uji is smaller 
than the radial frequency ujr by a factor A it is possible to 
fill the first N longitudinal states while the radial wave 
functions of the Fermions are still that of the ground state 
provided A^ < 1/A holds. 

From semi-classical theory or local density approxima- 
tion (LDA) it is well known (cf. e.g. g|), that a Fermi 
wave number kp = •\/(2Af — l)mujg/h can be associated 
with the one-dimensional Fermi gas of atomic mass m in 
a harmonic trap. It is noted that the condition N < 1/A 
is roughly in line with the standard estimate kp < l/l^ 
for a Fermi system which is confined to a transverse width 
It to be quasi one-dimensional provided the length It is 
identified with the extension Ir ~ ^Jh/{mLOr) of the ra- 
dial ground state wave function in the trap. 

Non-interacting Fermions in anisotropic harmonic 
traps have been studied recently p-[lC|] using exact and 



semi-classical methods. The thermodynamics of har- 
monically confined spin-polarized Fermions in any spa- 
tial dimension including a harmonic two-particle interac- 
tion has been studied in [ pd[ using the general approach 
p2| . The latter results are not available in closed form 
and require numerical evaluation. A finite series repre- 
sentation for the free energy of one-dimensional non- 
interacting spin-polarized Fermions confined by a har- 
monic potential has been given in |13| . 

In view of the feasibility to realize one-dimensional 
Fermions at ultra-low temperatures it seems worthwhile 
to supplement these works by studying the strictly one- 
dimensional case of non-interacting Fermions at zero 
temperature when a number of exact explicit results can 
be obtained. 

Interactions between spin-polarized identical Fermions 
are weak because the Pauli principle forbids s-wave scat- 
tering. On the other hand the theory of Luttinger liquids 
(cf. e.g., |1J] for a review) shows that even small inter- 
actions change a one-dimensional Fermi system substan- 
tially. Nevertheless, it is useful to have results for the 
non-interacting case to compare the effect of interactions 
with them. The results that we present below show fea- 
tures specific for one spatial dimension. 

In existing micro traps magnetic gradients of up to 
30 T/cm has already been realized [^ resulting in a pe- 
riodic motion of the trapped atoms on a time scale of 
micro seconds. Novel versions of micro traps based on 
micro fabricated current conductors achieve even higher 
gradients with an expected radial atomic oscillation fre- 
quency of above 1 MHz M . For the longitudinal oscilla- 
tion frequency 1 Hz appears to be a reasonable lower limit 
because time scales longer than a second gives rise to ex- 
perimental difficulties due to seismic and acoustic noise. 



Thus, the maximum value which is currently feasible for 
A is 10~^ and would limit the number of atoms inside 
the trap to about one million. The main experimental 
difhculty, however, is to fill the f 0^ states of the micro 
trap with a substantial number of atoms. Starting from 
an optically cooled sample of atoms with a phase space 
density of typically 10^^ Q a phase space compression 
of six orders of magnitude is required to completely fill up 
the wave guide. Such compression is possible with state 
of the art techniques of evaporative cooling |16|. Thus, a 
conservative estimation for realistic experimental condi- 
tions would assume a one component fully spin polarized 
Fermi-gas with a radial frequency inside the micro trap 
of 10^ Hz. The longitudinal frequency can be set at 10 Hz 
giving A = 10~^. Thus N = 10^ quasi one-dimensional 
Fermions can be accommodated inside the trap. Assum- 



ing Li atoms (in the hyperfine state 



TOc 



l/2,m, = 1)) 



the inverse harmonic oscillator length a according to 



— y'muji/h 



(1) 



is estimated as a « 8 • 10^ cm""'^ leading to a Fermi 
wave number kp « 10^ cm~^. 

Obviously the quasi one-dimensional Fermi energy ep, 
i.e., the energy of the highest occupied state without the 
radial contribution is 



ep = huji{N - -). 



(2) 



Under the above assumptions ep corresponds to about 
5 /iK and this temperature must be larger than the phys- 
ical temperature in order to achieve degeneracy of the 
Fermi gas. 

Another relevant quantity is the spatial extension of 
the inhomogeneous Fermi gas. The appropriate measure 
is twice that later given in equation (Q) and leads to 
a characteristic extension of 0.4 cm and to an average 
Fermion density of about 3 • 10** atoms per cm. The 
radial width 2/^ is about 3 • 10"^ cm. Thus the tonks 
gas limit |17| is avoided and the fermionic atoms can be 
treated as point particles. 

The exact quantum mechanical results usually give 
only small corrections to the corresponding LDA predic- 
tions. Some of them are, however, of qualitative nature 
and worth to point out. Among them are diverging den- 
sity oscillations near the classical boundary of the trap for 
large Fermion numbers and the general feature of Friedel 
oscillations [|l8[ of the density. 

The paper is organized as follows. Sec. II presents 
the basic theory. Sec. Ill discusses the relevant lengths 
and energy scales of the one-dimensional Fermi gas in 
the harmonic trap. In Sec. IV we compile the results 
for the zero temperature one-particle density distribu- 
tion. Sec. V is concerned with the Fourier transform of 
the density distribution. Sec. VI discusses momentum 
distributions and in the final Sec. VII we calculate the 
expansion of the particle density distribution when the 
trap is suddenly switched off. An Appendix summarizes 
mathematical formulae used in our calculations. 



II. BASIC THEORY 

We consider a gas of spinless non-interacting Fermions 
in one spatial dimension and trapped in a harmonic po- 
tential 



Viz) = 2^^^z^- 



(3) 



The Hamiltonian in second quantization and for the 
grand canonical ensemble is 



Ho = ^{hun -Ai)Crt< 



(4) 



n=0 



with one-particle energies fkUn — hLU({n + 1/2), n — 
0,1,.... The chemical potential is denoted fi. The 
Fermion creation operators c'^ and destruction operators 
c obey the fermionic algebra c„iC+ -I- c^Cm — 5in,n- This 
ensures that each (non-degenerate) energy level e„ = hujn 
with (real) single particle wave function 



^„W = y^^^e-"^^V2^,^(,,) 



(5) 



(normalized according to {m\n) = 5m,n) is at most 
singly occupied. The intrinsic length scale of the system 
is the oscillator length I — a~^ where a is defined by (H). 
Hn denotes a Hermite polynomial. 

We consider the spatial density of one-dimensional 
Fermions in the harmonic trap, i.e., the one-particle dis- 
tribution function 



n{z;T,fj.) = Tr pip^{z)'ilj{z). 



(6) 



In (pp the operator tp(z) destroys a Fermion at position 
z. It can be expanded as 



1p{z) = y^ V'n(z)Cn. 



The density operator is 



(7) 



(8) 



with Z = Tr exp[— /3_ffo]- A standard textbook exer- 
cise then gives 



no 



{z-T,n) ^ Y^ i>l^{z)pm{T,^.) 



m=0 



where 



p™(r,Ai) = {e'^(^"- -'') + !} ' 



(9) 



(10) 



is the thermal occupation number of the single particle 
state ipm- 

The present paper deals with the case T ^ Q when a 
number of analytical results are available. The impor- 
tant simplification results from the fact that for T ^ 
the first TV levels are completely filled while all others are 
empty, i.e., 



Pm{T^O,^i)^eiN-l-m) 



and /i becomes the Fermi energy ep'- 



H^ ep = fuji{N - -) 



(11) 



(12) 



The density no{z;T —* 0,fj,) which we henceforth de- 
note uq^z, N) takes on the form 



N 

o 
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FIG. 2. Particle density distribution functions in units of 
the inverse oscillator length a for N — 5 and A*' = 6 Fermions 
in a one-dimensional harmonic trap and at zero temperature. 
The added Fermion resides in the area between the two curves. 
The density oscillations near the center can be identified as 
Friedel oscillations. 



7V-1 



no{z,N)^ 51 ^"(^ 



(13) 



n=0 



The zero temperature case is depicted in Fig. 1. 
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FIG. 1. A'' = 10 non-interacting spinless Fermions filling 
the lowest single particle levels in the harmonic trap at zero 
temperature, cf denotes the Fermi energy and Lp the half 
width of the Fermi system. 



Equation ( |13D is the main object of the present study. 
Under the condition fc^T <C cp it correctly describes the 
density of non-interacting Fermions in a harmonic trap. 

Fig. 2 shows the density profile with the characteristic 
ripples on top. 

This is in contrast to an infinite Fermi gas (or one with 
periodic boundary conditions) where the density is homo- 
geneous. The ripples appear here as a finite size effect. 
In the center of the trap they will be identified below as 
Friedel oscillations [fsl. 



III. LENGTHS AND ENERGY SCALES 

In this Section wc summarize the relevant scales of 
the one-dimensional Fermion gas in the harmonic trap. 
They are expressed in terms of the basic quantities m, 
uji, N. One of them is clearly the Fermi energy ep ac- 
cording to (p^ ) . At the Fermi energy the filled Fermi sea 
has a spatial extension 2Lp according to mujfL'jp/2 = 
huJi{N- 1/2) or 



= -V2iV - 1 = L„= 



Af-l, 



(14) 



a quantity frequently appearing later. The positions 
z = ±Lp are classical turning points for a Fermion with 
energy ep. The length Lp is the largest length of the 
problem followed by a~^ which is associated with the 
zero point en ergy . 



Equation (a.8) of the Appendix shows that a wave 
function ipn behaves as a standing wave with wave vec- 
tor kn — ay/2n + 1 in the middle of the trap provided 
n 3> 1. At the Fermi energy the wave number be- 
comes kp = a\/2N — 1. Together with ( |l^ ) this leads 
to ep — h'^kp/2m as suggested in g] for the three- 
dimensional anisotropic case. The Fermi wave number 
kp^ is the shortest length scale of the problem and the 
Fermi energy the largest energy. 

What is the relation between kp and the particle den- 
sity? In a one-dimensional Fermion gas of spatial exten- 
sion 2Lp with periodic or open (infinite potential well) 
boundary conditions the relation is in both cases 



k]p' = TTOo 



2I7f 



(15) 



In the present inhomogeneous situation kp increases 
as iV^/^ because the width 2Lp of the trap also increases 



as N'^/'^. However, we can discuss peak density ng and 
average density uq (or even higher moments of no(z, N)). 
The peak density i s cl early found near z = 0. Using the 
asymptotic result (|a.9[) 



{l-{-iy^ cos2kFz) (16) 



fcp 1 



gives 



,(p) 



kp 

n 



(17) 



The sign ~ denotes here and further on an asymptotic 
correspondence for N :$> 1. Note that this asymptotic 
limit does not imply the semi-classical approximation. 

From ( [l7| ) it is seen that the usual relation (|l^) be- 
tween kp and the one-dimensional particle density refers 
here to its peak value near the center. 

It is more difficult to discuss the average density hq 
since an averaging length is needed. Equations ([2^) or 
(pl|) only give the obvious sum rule 



dzno{z,N) = N. 



(18) 



We thus resort to the semi-classical approximation 
where the local density is given [|| by 



nsc{z,N) 



_ kp{z) 



k 



^-/-fe) 



(19) 



which is zero outside \z\ < Lp. The corresponding av- 
erage density clearly is uq = kp/A and is only slightly 
smaller than the peak density. 

The sum rule for Usc gives N — 1/2, i.e., half a Fermion 
is missing under the curve Ugc- This is due to the ne- 
glect of Fermion density leaking out of the classical re- 
gion \z\ < Lp by tunneling. One might conclude that the 
number of Fermions in the oscillations is about one half. 
This is not correct for large Fermion number TV when 
the difference between no{z,N) and nsc{z,N) near the 
boundaries becomes significant due to increasing oscilla- 
tions in the exact density as detailed in the next Section. 

Finally, using formula (a.9) immediately allows the 



identification of the ripples in nQ{z,N) near the center 
with the well known Friedel oscillations pq ] of wave num- 
ber 2kp around an impurity in the degenerate Fermi sea. 
In a naive interpretation these oscillations result from the 
superposition of incoming and reflected parts of the up- 
permost wave function which both have a wave number 
kp near the center of the trap. A more subtle interpre- 
tation refers to the inherent instability of the degenerate 
free Fermi gas towards static longitudinal perturbations 
of wave number q = 2kp. A well studied example are free 
electrons (cf. e.g., [M). While in three dimensions only 
a logarithmic singularity in the derivative with respect 
to q appears in the susceptibility it becomes a logarith- 
mic singularity in one dimension due to perfect nesting. 



This causes charge and spin density instabilities when 
backscattering interactions are present. In bounded Lut- 
tinger liquids (cf. ||2^,|2l| ) the interactions modify the di- 
vergence of the density oscillation near the boundary. 

But even without interactions breaking of transla- 
tional invariance by inhomogeneities like impurities and 
boundaries trigger density oscillations of wave vector 
2kp. In one dimension this effect is well known for non- 
interacting Fermions with open boundary conditions. 
While in one dimension the effect is most pronounced it 
is, nevertheless, possible to identify the oscillations in the 
isotropic density calculated in |^ as three-dimensional 
Friedel oscillations. 

Since the one-dimensional Friedel-oscillations contain 
only about one atom it will be difficult to detect this ef- 
fect experimentally. However, it is conceivable to use an 
array of shorter micro traps each filled with a reduced 
number of atoms. The oscillations within each trap then 
add up and lead to a total effect that is enhanced by 
the number of traps. Using micro fabrication techniques 
it should be possible to combine 100 traps on one sub- 
strate leading to a signal that may become within reach 
of advanced imaging techniques. 

Friedel oscillations though difficult to observe are a 
fundamental property of the degenerate Fermi gas which 
eludes the semi-classical approximation. 



IV. ONE-PARTICLE DENSITY DISTRIBUTION 



With the help of (a.3) the summation in dl3) can be 
performed for any z and N with the result: 



noiz, N) = N^%{z) - y/N{N + l)V'jv+i(^)V'JV-i(z). 

(20) 

Using the recurrence relations for the wave functions 
'0n (cf. Appendix) this expression can be brought into 
another useful form 



no 



(z,iV) = iVV4-i(-z) - VN{N-1)iPn{z)4>n-2{z). 



(21) 



This shows that the density distribution must be a 
polynomial of order iV — 1 in a^z^ times the exponential 
exp{— a^z^} since the density is an even function of z. 

The formulae (g^) and ( pi| ) admit a number of exact 
conclusions as well as some remarkable asymptotic results 
with respect to the Fermion number N. 

Differentiating ( po|) and ( pl|) w ith respect to z and us- 
ing the recurrence relations (a.l) and (|a.2|) gives 



dnoiz,N) 
dz 

d^no{z,N) 
dz"^ 



-aV2NipN{z) il^N-i{z), 



2a^N[ij%iz)-ib%_^{z)]. 



(22) 



This shows that the density distribution no{z, N) has 



• A'' maxima at the N zeros zl ' of tp^ {z) , 
(z.= l,...,iV), 

• A^— 1 minima at the iV— 1 zeros zl oiipi\f-i{z), 
(i. = l,...,7V-l). 

As a consequence the minima of no{z, iV+ 1) touch the 

maxima of no{z, N) at the points zl . This is shown in 
Fig. 2 for iV = 5. The area between no{z, 6) and no{z, 5) 
contains precisely one Fermion. In this way the Pauh 
exclusion principle is optimally implemented. 

The above considerations also show that about half a 
Fermion is contained in the ripples of the density distri- 
bution. The density at the maxima is given by 



AN) 



(^')^ 



and at the minima it is: 



(23) 



(24) 



Due to the knot theorems |p2| the topological features 
inherent in the above statements carry over to arbitrary 
concave potentials. Thus counting the number of maxima 
of the density distribution gives the number of Fermions 
in any concave trap. 

We now come to asymptotic results for iV ^ 1. In 
practice, A^ « 20 is a good lower bou nd. In the asymp- 
totic region the powerful formu la ( a.4) is available for the 
full range \z\ < Lp- Inserting ^A) into ( pl| ) gives 

no(z, N)^kp{{l + ^) (zl^Ai^(i^_0 (25) 
- (1 + ik) r^'T'^'^T Ai(i^)Ai(t^_2)l . 



The functions t„(z) and (t>n{z) are defined in (a. 5) and 
( |a.6D . The advantage of this formula lies in the fact that 
the indices of the wave functions moved into the argu- 
ments of the Airy functions. 

Evidently, the positions zl of the maxima are now 
found from 



Ai 



i(^A'(4^^) 



0. 



(26) 



We are interested in the positions of the last maximum, 
i.e., those lying in the neighbourhood oi Lp. The asymp- 
totic expansion of t^ in the region z < Lp ~ (2A^)^/^ is 



tN ~2iV2/^ ( 1- T— 
Lp 



One also finds 



t{N-l)±l ^ tN-lT N ^'^. 



(27) 



(28) 



Provided N-^'^ is much larger than unity this leads to 
(the prime means the derivative) 



no(zW,iV) ^ aV2N'/'Ai' (iA.-i(zW)) (29) 

~aV2 7Vi/6Ai'(tA.(zW))2. 

Specifically, we consider the last maximum at z]y . It 
corresponds to the first zero of the Airy function Ai(t) 
which is at tAr(z]y ') = —2, 33... . 

This gives 



Zm ' ^ Lp { I- TTTTTT ] ^ Lp - Axn- 



-N 



Note that 



Ar2/3 



1.17 ^ 1.65 



iV2/3 " aN^^^' 
while the density at the maximum is 

no(z^\7V)~0.7a7Vi/^ 



(30) 



(31) 



(32) 



In the same way the distance I = Zj^ — zj^_^ between 
the last two maxima can be calculated. It determines the 
smallest local wave number fc(™™) = 271/1 which is found 
to be 



fc('^ 



2Tra0.8N^^^. 



(33) 



It needs more than A^ = 50 Fermions to make fc*^™*"^ 
less than half the maximal wave number fc('"°=^) = 2kp 
appropriate for the central part of the trap. 

When we define a shrinking region S{N) near Lp ac- 
cording to 



Aa; = Lp — z ~ 



f 



with - say - / varying from 1 to 10 then we have 



(34) 



(35) 



which is independent of N . 

In S{N) the density ( |25| ) can be drastically simplified 
to read 

no(z, N) - aV2N^'^ lAi{~tNMz)f (36) 

-Ai(iAr_i(z))Ai"(tA,_i(z))}. 

Thus (^9|) leads to a self similarity of the tails of the 
density for z G S{N) defined above: Calhng no(Aa;, A^) — 
no{z = Lp - Ax,N) the graph of no(Aa; e S{N2),N2) 
maps precisely onto that of no{Ax S S{Ni),Ni) when 
the density is rescaled according to {Ni/N2y^^ and 
the position Ax according to (A^2/Afi)^'^- Note that 
nsc{Ax,N) also satisfies the scaling in S{N). 



Finally we exploit the approximation (a. 7) which holds 
for large values of (-t), i.e., slightly away from the classi- 
cal turning points z = :kLp. We are interested in results 
for \z\ < Lp which better the result (a.£) valid in the 
very middle of the trap. 

Using the expansion 



(N-l)±l 



yjv-1 



± 



2N- 



1 VL 



(37) 



which requires distances Az/a away from the bound- 
ary to be much larger than N~'^'^ we find 



jio(z,N) = n,c{z,N) + 



2-,tLf^/1-(z/LfY 



(38) 



2iTLf 



\U2N-1)(^1-{z/Lf)^ 2/LF-arccos(2/LF)')| 



Equation ( |38| ) separates the slowly varying background 
Use from an increasing and spatially oscillating part due 
to quantum effects. 

It is also seen that the oscillating part increases to- 
wards the boundaries z = ±Lp. A naive extrapolation 
would give an envelope 



E{z) 



1 



TTLpy/l - {Z/Lpf 



(39) 



of these oscillations which formally diverges near Lp 
as (1 — [z/LpY)~^/'^ . In view of the range of validity 
of (|3^) this is, however, unwarranted. Nevertheless, it 
raises the question how the oscillating part diverges at 
the boundaries when N diverges. 

N non-interacting Fermions in a one-dimensional box 
of width L confined between z = and z = L^ and 
with infinite barriers (open boundary conditions) have 
the density distribution (for A^ ^ 1) 



no(z,A) 



^(0) 



1 — cot 



7rz\ sin2fc 



L 



(0). 



2N 



(40) 



Using (|l|), i.e., kp 
z <^ L according to 



(0) 



■kN/L this gives an envelope for 



E{z)^-^—^. 



(41) 



We conjecture that in our case of a soft boundary the 
limiting behaviour near the right boundary for very large 

N is 



E{z < Lp) 



1 



2TTLp{i-z/Lpy(^y 



(42) 



There is numerical evidence for 5{N ^ oo) ^ 1, how- 
ever, in a very slow approach {S{N) « 1 — 1/ln A^). This 
would imply that the integrated absolute fluctuations 



l-LF 

5N ^ dz\no{z,N)-nsc{z,N)l 

J-LF 

i.e., the number 5N of Fermions in the density oscilla- 
tions diverges logarithmically with N as it does for ideal 
Fermions in a box. The mathematical problem in clar- 
ifying this point lies in the enormous difficulty in sub- 
tracting out the oscillating part in the boundary region 
z -^ Lp{N ^ 1) which is outside the approximation 

@- 

In case of interacting spinless one-dimensional 
Fermions in a box it is known that the exponent 6 is 
given by the coupling constant K pfl] . 



V. DENSITY PROFILE IN FOURIER SPACE 

With possible application to optical detection we dis- 
cuss the Fourier transform of no{z, N) 



Fno{k,N) 



dze'^^'noiz^N). 



(43) 



It can be evaluated exactly in the following way: The 
integral in [E3[ can be converted to the form 

jT^dze^'^^i^U^.)M^') (44) 

— c^t4- f -k^ \(n-m)/2 fm\ t (n-rn) , k'^ s 

(n > m) where -L„ denotes a Laguerre polynomial. 
Applying (B3) to ( pT| ) and using recursion relations for 
Laguerre polynomials |pl| gives 



Fn^ik.N) 



_J^ (1) , K . 



(45) 



The Fourier transform of the semi-classical expression 
(ES) can also be given in closed form involving a Bessel 
function: 



Fn,,ik,N)='^Jiij-^). 



kp 
~k 



(46) 



Note the sum rules 

dkFno{k,N) = 2kp 

and the limits 



dkFn,,{k,N), (47) 



i^no(A:^0,A^) = A^, Fn^dk -^ 0,N) ^ N - -. (48) 

The basic difference between the exact result (HS) 
which takes care of the ripples in the density profile and 
the semi-classical form (Hfl) is a hump somewhat below 
the wave number 2kp as shown in Fig. 3. 




10 



15 
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FIG. 3. Part of the Fourier transformed density distribu- 
tion function for A'^ = 20 non-interacting Fermions at zero 
temperature. Note the small hump near the wave number 2kF 
due to Friedel oscillations which provide the shortest length 
scale of the problem. Dotted line is the semi-classical approx- 
imation lacking that feature. 

For larger wave numbers Fno{k, N) drops to zero while 
Fnsc{k, N) shows a multitude of oscillations similar to 
that produced by a slit of spatial width 2Lf- For 
k ^ 2kp, however, the exact result and the semi-classical 
approximation agree very well. 



VI. ONE-PARTICLE MOMENTUM 
DISTRIBUTIONS 

Even for a confined system one can define a momentum 
density distribution by 



p{k) - /^+ 



Cfc Ck) 



(49) 



The operator Ck annihilates a Fermion with (contin- 
uous) momentum hk. It can be decomposed into the 
fermionic annihilation operators for the harmonic oscil- 
lator according to 



n=0 

with the transformation function 

i™ k 

ft = — V'm(z = -^). 



(50) 



(51) 



The momentum density po of non-interacting 
Fermions in the harmonic trap and at zero temperatures 
thus becomes 



Po(fc; TV, T = 0) = E (f';Xft Q{N~\- m) (52) 

N-1 



leading to the remarkable result 



PO 



(fc;iv,r = o) = ^E^™(^ = ^/"') 



(53) 



m— 



= — no(z = k/a^; N,T ^ 0). 

The momentum density is isomorphic to the particle 
density with kp replacing Lp. Obviously it satisfies the 
general sum rule 



f 

J — ( 



dkp{k;N,T)^N. (54) 

Alternatively, we can study the momentum probability 
P{k;N,T^O)= / dze''^^(Vi+(z)V;(0))W. (55) 



which is also appropriate for symmetric confining po- 
tentials centered at z — 0: 

For non-interacting Fermio ns a t zero temperature in 
the harmonic trap we can use (a. 3) to find: 



(V'+(z)V(0))r^ = ^e 



" ~a^z^/2 



(56) 



El 

M 



dN,2M+iL'^^i\a^z^) + SN^iL^jjliia^z^)} 



In the limit N ^ 1 this reduces to the simple expres- 
sion 

(^+(.)^(0))r'-^+o(-^). (57) 



The corresponding discrete momentum distribution for 
the harmonic oscillator is the well known step function 



Po{k; N:$>1,T^0)^ e{kp - k) 



(58) 



with k = kn = ay/2n + 1. 

Some general remarks may be useful: The centered 
momentum distribution (55) can be expressed as 



P{k-N,T) 



dk' {clck')- 



(59) 



In case of translational invariance the integrand in (p9) 
becomes 



{ctck') =P{k-N,T)5(k-k'). 



(60) 



m=0 



Under periodic boundary conditions both definitions 
(|i) and (11) coincide. 



VII. FREE EXPANSION OF PARTICLE DENSITY 

Expansion of a particle cloud is an important tool to 
investigate Bose-Einstein condensates (c.f. 24, y]). De- 
tailed theories are available for this expansion based on 
the Gross-Pitaevskii equation p4|-p8[. In the simplest 
case the trap is suddenly switched off and the conden- 
sate expands freely. Non-interacting Boson condensates 
display universal length scaling [g5[ in all spatial dimen- 
sions. The same has also been shown recently for non- 
interacting Fermions [|lO| . 

When particle interactions dominate the kinetic energy 
the transverse scaling function in a highly anisotropic 
trap is 



br{t) = y/T+Z^, 



(61) 



while the longitudinal expansion is more complicated 
[26j. This was confirmed experimentally in full detail in 
|29|.|0|. 

In accord with llfl] we find that a freely expand- 
ing degenerate one-dimensional gas of non-interacting 
Fermions behaves according to mWj with ue replacing Ur- 
The calculation is fully quantum mechanical and supple- 
ments the approach in [^ . 

The quantity to be calculated is 



n{z,t) ^TTp{t)^ij-^{z)iP{z). 



(62) 



where p(t) is the density operator of the freely expand- 
ing gas. It is given in terms of the statistical operator p{0) 
immediately before the trap is opened at time t = by 



m 



— e^*"^""*/''/5(0) e^^°°^/^ 



(63) 



The free expansion of non-interacting one-dimensional 
Fermions is governed by the Hamiltonian 



Hn 






dk — c,Ck. 



(64) 



The operators c^ and c^ were introduced in conjunc- 
tion with (E9). Equation (p3) can also be written as 



with 



n(z, t) = Tr /5(0)V'+ (z, t)iP{z, t) 



^(z,t) = e'"««^'^{z)e-'"°°^'^ 



^j^^r(kz-uut) 



(65) 
(66) 



and 



We now use (| 



/QjT 



Wfc 



Cfe 



hk"^ 



2m 
and (P) to find 



(67) 



no[ 






27r 



E (~l)'"+"(/™)*/n'Tr/5(0> 



n+. 



m,n^O 



For a harmonic trap which is initially in thermal equi- 
librium the statistical operator /5(0) is (@) and at zero 
temperature 



Tr pc+c„ = dm,n 6 (iV - 1 - to) 
holds. We thus find 

1 



no{z,t;T ^ 0) 



27ra^ 



dkdk' 



(69) 



(70) 



[m=0 ^ a j 

The summation in curly brackets can be performed 
using (a.3). In order to proceed it is convenient to go 
over to the Fourier transform and write out the oscilla- 
tor eigenfunctions in terms of Hcrmitc polynomials. This 
leads to 



FnQ{ki,t;T = 0)= dz e"^^'' m{z,t-T = Q) (71) 



TT 2Nkl{N ~l)\ 



e 2 



'\{l-iujet) 



dk' e- 



'-^ ^(1— 2a;^t) 



U^fl±^) HN-ii-) -{N^{N- 1)) 
a a 



Using [ p2| the integration can be performed giving 

Fnoik, t; T = 0) = e-l&(i+-^^*^) L^nU^^ + u^^))- 

(72) 

This formula is isomorphic to ( [iq ) with the inverse 
length a being replaced by the rescaled value 



a^{l+ uj'jt^)-^/'^a = a/b{t). 



(73) 



Since Fno(k,t) and 710(2:, t) are related via a Fourier 
transformation the final result is 



no(.,t;T.O)^^no(^,iV). 



(74) 



Thus free longitudinal expansion proceeds via a simple 
length rescaling involving the factor b{t). In the course 



of time the initial density distribution Fig. 2 decreases 
and broadens according to the factor b{t) but preserves 
its topology including the Friedel oscillations which cor- 
respondingly increase their wave length. 

In this picture it is assumed that the Fcrmions remain 
one-dimensional during the expansion. If the transverse 
confining fields are also removed transverse expansion in 
any of the two equivalent transverse directions will also 
proceed according to (174) taken for A^ = 1 and with the 
scaling function (|6l| ) . This follows simply from the obser- 
vation that each of the two ground state wave functions 
for the transverse directions correspond to a single one- 
dimcnsionally confined Fermion with ujr in place oi ug. 



VIII. SUMMARY 



■lpn+l{zi)-4'niz2) " '4'nizi)lpn+liz2) 



a{zi - Z2) 



which is a conversion of a formula in ]33[ |. In (a. 3) the 
limit zi — > Z2 can be performed and the resulting deriva- 
tives c onv erted into harmonic oscillator wave function 
using (a^) and (aJ). This leads to (PQ). An alternative 
derivation applies induction to ( |2l| ) which is o bvio usly 
true for A^ = 1 utilizing the recurrence relation (a.l). 

A very useful asymptotic (n 3> 1) expression for the 
wave functions can be extracted from |31| (Chap. 19.7): 



i^n{z) 



2 

'a I - 

n 



with 



1/4 



(-in) 



1/4 



sm 



1/2, 



-Ai(i„) 



(a.4) 



We have calculated exactly one-particle properties of 
non-interacting one-dimensional Fermions in a harmonic 
trap. These are the particle density distribution includ- 
ing its free expansion when the trap is switched off and 
also two momentum distribution functions. The exact 
calculability can be traced back to two specific mathe- 
matical features of the eigenfunctions of the harmonic 
oscillator namely that finite sums of bilinear expressions 
can be performed and that Fourier transformation essen- 
tially reproduces an eigenfunction. Friedel oscillations in 
the particle density and its analogue in the momentum 
distribution as well as diverging density oscillations near 
the classical boundary are basic features of the degener- 
ate one-dimensional ideal Fermi gas. 
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IX. APPENDIX 



t-n 



3 f n 1 , ,„ , . „ , X 

2 ( 2 +i)(2</'n-sin20„) 



2/3 



and 



C0S( 



z 



(a.5) 



(a.6) 



Ai is the Airy-function which oscillates for negative ar- 
guments. There is a continuation to positive arguments 
(the tunneling region) which we will not discuss. 

The tilde ~ denotes asymptotic expansion for large n 
including all prefactors. Inside the trap, i.e., away from 
the classical borders, the form 



1pn{z) 



a I 2 

TT V 71 



1/4 



1 



sm 



1/2, 



(a.7) 



Tl\\ TT 

*^°^ *! I o + 7 (sin2(/)„ - 20„) + - 



The Appendix is a compilation of some mathematical 
formulae used in the derivation of the results given in 
the main part of the paper. An important role is played 
by the recurrence relations for Hermite polynomials (cf. 
e.g., pHI). Here they are given as recurrence relations for 
the complete harmonic oscillator wave functions ipn{z)- 
These are 



is useful. It results by means of the asymptotic ex- 
pansion of the Airy function Ai(t) for —t^l. In 
the limit \z\ <C Ln a further simplification occurs since 
(j)n — > Ti'/S — z/Ln- This leads to (note: n ^ 1) 



ipn 



1/4 



~~2 ) '^OM 



n^ri.Z - 






(a.8) 



Vn+ lijjn+i{z) - azV2^n{z) + ^/nipn-iiz) = 0, (a.l) 



which is used in Sec. IV. The corresponding Fermion 
density well inside the trap is 



-—ip„{z) + a^zipniz) - aV2n'ipn-iiz) = 0. (a.2) 
az 



no(z,7V) 



1 



2ttLf 



(1 - (-1)'''cos2A:f^ 



(a.9) 



The summation of the finite series ( |13[ ) is accomplished 
by means of 



y^ ■4^,n{zi)ip„i{z2) 



(a.3) 



m=0 



Here a small systematic error l/(2TrLp) of this approx- 
imation for \z\ ^ Lp has been subtracted to bring ( |a.9| ) 
in line with the exact result. 
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